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Abstract. In this paper we proposed a verified numerical method for deriving a sharp inclusion 
of the Sobolev embedding constant from 77,) (P) to L P (P) on bounded convex domain in M 2 . 
We estimated the embedding constant by computing the corresponding extremal function using 
verified numerical computation. Some concrete numerical inclusions of the constant on a square 
domain were presented. 
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1 Introduction 

We are concerned with the best constant C p (P) in the Sobolev type inequality satisfying 

IMIlp(Q) — C'p(P) Htfi(n), Vu € Hq (P), (1) 

where P C M n (n = 2, 3, ■ ■ •) and 1 < p < oo if n = 2, 1 < p < (n + 2)/(n — 2) if n > 3. 
Since Sobolev type inequalities are important in studies on partial differential equations, there 
have been a lot of works on such inequalities and their applications, e.g., [1, 2, 3, 4, 5, 10, 13, 
18, 21, 20, 23]. The classical Sobolev embedding theorem has been well known. Moreover, a 
formula giving the best constant in the classical Sobolev inequality on R n was independently 
shown by Aubin [1] and Talenti [23] in 1976 (see Theorem A.l). Since all elements in H ( \ (P) 
can be regarded as those in H l ( P) by zero extension outside P, we can obtain a rough upper 
bound of C p ( P) for a general domain P C M n using the formula (see Corollary A.l). One can 
find another estimation formula in [18] (see Theorem A.2). 

In this paper we will propose a numerical method for deriving a verified sharp inclusion of 
the best constant C p (P) satisfying (1) for bounded convex domain PcR 2 , e.g., we proved the 
following proposition by our method: 

Proposition 1.1. The smallest values of C p (P) (p = 3,4,5) satisfying (1) are enclosed as fol¬ 
lows: 

C 3 (P) € [0.25712475017617, 0.25712766496560]; (2) 

C 4 (P) € [0.28524446071925, 0.28524446071939]; (3) 

C 5 (P) € [0.31058015094169, 0.31067136032829], (4) 

where P = (0, l) 2 . 

Hereafter, we replace the notation C p (P) with C p +\ (P) for notational convenience. The 
smallest value of C p + 1 (P) can be written by 

Cp+ 1 (P) = sup 4> (u) , (5) 
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where 4> (u) = IM|/>+i(q) / IMIffi(n)- It is well known that C p+ \ (12) in (5) is finite and realized 
by an extremal function u* € Hq (12) (see, e.g., [6]). The critical point problem for <f> is reduced 
to finding weak solutions to the following problem: 


—Au = u p 
u > 0 
u = 0 


in 12, 
in 12, 
on <912 


( 6 ) 


This problem has only one solution if 12 C K 2 is bounded and convex [11]; therefore, in this case, 
we can derive an inclusion of C p+ \ (17) = ||'U*|| L p+i( 5 - 2 ) / ||^*||#i(n) by computing the solution to 
(6) with verification. 

There is a number of numerical methods for verifying solution to semilinear elliptic boundary 
value problems (e.g., in [14, 17, 18, 22]) and related works, e.g., [15, 24], These methods enable 
us to obtain a concrete ball in the senses of ||V-|| l2 (q) and ||‘llz,°o(0) containing exact solution 
and can be applied to the problem: 

J — An = |u| p_1 u in 17, , , 

1 u = 0 on <917. 


However, positiveness of the solution to (7) should be proven to derive a verified solution to (6). 
One of the main contributions of this paper is to propose how to verify the positiveness. 

This paper consists of the following sections: In Section 2 and 3, we propose a method for 
proving positiveness of a solution to (7) and a method for estimating the embedding constant 
C p +1 (17), respectively. In Section 4, some numerical examples are also presented, which lead 
Proposition 1.1. 


2 Verification method for positiveness 

In this section, we propose a sufficient condition for positiveness of the solution to (7), which will 
be summarized in Theorem 2.1. This theorem enables us to numerically check the positiveness. 
Let us first introduce the following notation to be used throughout this paper: 

• define N := {1, 2,3, • • • }; 

• let B ( x , r ; || • ||) be an open ball whose center is x and whose radius is r > 0 in the sense 
of the norm || • ||; 

• denote its closure by B (x,r ; || • ||); 

• let L°° (17) be the functional space of Lebesgue integrable functions over 17, s.t., \u (.x)| < 
oo (a.e. x € 17) with the norm ||u|| L oo(^- ) := esssup{|u {x)\ \x € 17}; 

• let H 1 (17) be the first order L 2 Sobolev space on 17; 

• let Hq (17) := jtt G H 1 (17) : u = 0 on dQ in the trace sense} which endowed with the 

norm IHIffi(n) : = II v 'IIl 2 (Q)- 

Let us remark that, for Lebesgue integrable functions, we omit the expression “almost every¬ 
where" for simplicity, e.g., we denote u > 0 in the place of u(x) > 0 a.e. x € 17 for u € Hq (17). 
Let us introduce the following lemma for proving Theorem 2.1. 

Lemma 2.1. Let 17 be a bounded domain in M n . All weak solutions u € Hq (17) to (6) satisfies 
esssupju (x) p_1 | x € 17} > Ai, where \\ > 0 is the first eigenvalue of the following problem: 

(V“, Vu) L 2 ( q) = A (u, u) L 2 ( q) , Vu € Hq (17). (8) 
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Proof. Let 0i > 0 (0i ^ 0) be the first eigenfunction corresponding to Ai, which satisfies 


We have 


/ u p (x ) 4>i (x) dx = Ai u ( x ) 0i ( x) dx. 

Jci Jn 


/ u p (x) 0i (x) dx 

Jn 


Jn 


{u (x) 0i (x)} dx 


< M u / u (x) (x) dx 

Jo. 

= A f l M u f u p (x) 0i (x) dx, 
Jn 


where M u := esssup{u (x) p 1 | x € 0}. Positiveness of / u p (x) 0i (x) dx implies M u > Ai. □ 

dn 

Using lemma 2.1, we are able to prove the following theorem, which gives a sufficient condition 
for positiveness of the solution to (7). 

Theorem 2.1. Let C M n be a bounded connected domain. If a solution u € C 2 (LI) to (7) 
is positive in a subdomain II C 0 and if sup{u_ (x)^ 1 | x € f2} < Ai (LI), then u > 0 in the 
original domain 17, where \\ (Ll) > 0 is the first eigenvalue of the problem (8) and U- is defined 
by 


/ \ ( -u(x), u(x) < 0, 

U- (x) := 7 . n 

[0, u (x) > 0. 

Namely, u is also a weak solution to (6). 

Proof. Assume that u is not a positive solution in Ll, which ensures there exists a domain 
Ll" C Ll\Ll' such that u < 0 in Ll" and u = 0 on dLl"\ Hopf’s lemma makes it impossible for a 
subdomain LIq C Ll such that u = 0 in LIq to exist. Therefore, the restricted function v := —u 
can be regarded as a solution to 

( —Av = v p in Ll", 

< v > 0 in Ll", 

( v = 0 on dLl". 

From Lemma 2.1, we have 

sup it_ (x) p_1 > sup v (x) p ~ l 

xeVL x£tt" 

> Ai (17") . 

Since 17" C 17, we have Ai (17") > Ai (17). Hence, u is a solution to (6) if sup{u_ (x) p ~ l \ x € 
17} < Ai (17). □ 

Remark 2.1. i) Since Hopf's lemma requires the regularity u € C 2 (Ll), we also need this 
regularity to prove Theorem 2.1. A weak solution u € idg (17) to (7) is in C 2 (17), e.g., when 17 
is a bounded convex domain with piecewise C 2 boundary (see, e.g., [7]). 

ii) The first eigenvalue of the problem (8) can be numerically estimated by, e.g, the method 
in [12], which enables us to concretely evaluate the eigenvalues of (8) on polygonal domains in 
M 2 . 


3 



3 Estimation method for embedding constant 


In this section, we propose a method for estimating the embedding constant C v +\ (Q) defined in 
(5). The following theorem gives concrete estimation of the embedding constant from a verified 
solution to ( 6 ). 

Theorem 3.1. Let ft C M 2 be a bounded convex domain. If there exists the solution to (6) in 
a closed ball B(u,r ; || • ||#i(n)) with u € Hq (O) satisfying ||w||j/i(Q) > 2r, then the embedding 
constant C p+ \ (12) defined in (5) is estimated as 

IHIz>+i(f2) / ^ II^IIlp+1(Q) 

-p-f]- < Cp+1 (ft) < ]prii- 7T- 


Proof. It is obvious that |HIlp+i(Q) / IHIhI(q) is a lower bound of C p+ \ (12). The solution to ( 6 ) 
is unique when ft C R 2 is bounded and convex [11]; therefore, the ratio |M|£,p+ 1 (q) / |M|#i(Q) 
is maximized by the solution u to ( 6 ). Let us write the solution to ( 6 ) as u + rv with v € 
Hq [ft) , [Ml H^(n) = 1- Then, we have 


C p+ 1 (ft) 


u + rv\\ L p+i^ IHIlp+i(T2) + r Cp+i (^) 
ll“ + HI/ri(n) ~ ll“llH 0 i(n) ~ r 


that is, 


(lNlf/i(n) — 2r ) Cp+i (^) ^ H“llLP+ 1 (f2) • 


Hence, when ||&||#i(Q) > 2 r, ||h||£,p+i(Q) /(W^Wh 1 (n) ~ 2r ) becomes an upper bound of C p+ \ (ft). 

□ 


Remark 3.1. Theorem 3.1 can be naturally applied to the case that n > 3. For example, if ft 
is a convex symmetric domain in W 1 (n > 3) and 1 < p < (n + 2)/(n — 2), ( 6 ) has only one 
solution [ 8 , 16]. 


4 Numerical example 

In this section, we present some numerical examples of proving positiveness of a solution to (7) 
and estimating the corresponding embedding constant, which will lead Proposition 1.1. 

All computations were carried out on a computer with Intel Xeon E5-2687W 3.10 GHz, 512 
GB RAM, CentOS 6.3, and MATLAB 2012a. All rounding errors were strictly estimated using 
toolboxes for verified numerical computations: INTLAB version 6 [19] and KV library version 
0.4.16 [9]. Therefore, the accuracy of all results is mathematically guaranteed. 

We first treated the case that p = 3 and ft = (0, l) 2 , which corresponds to the critical point 
problem for the embedding constant C 4 (ft). In this case, (7) has an infinite number of solution, 
while it has only one positive solution [11], We computed an approximate solution to (7) with 
the Fourier basis <f>ij := a t j sin( 7 T'<x) sin(njy), 1 < i,j < N, atj € M. Moreover, we proved 
existence of the solution to (7) in a i^Q-ball and L°°-ball whose center is the approximation 
using the method in [17] with the method in [24] 1 . Figure 1 shows an approximate solution u 
to (7) such that existence of the positive solution to (7) in the a neighbourhood of u was proven. 
Table 1 shows verification results in this case: The first eigenvalue Ai (H) of ( 8 ) is 27r; therefore, 

^The method in [24] is for estimating the inverse norm of elliptic operators, which take a central role in verified 
numerical computation for partial differential equations. 
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sup{it_ (x ) p ^ 1 | x € fl} < Ai (0) holds, this means positiveness of the solution to (7) is proved, 
in all cases in the table. Moreover, the last row shows intervals containing C 4 (f2), e.g., 1 . 2345 ! 
represents the interval [1.23456,1.23789]; the case that N = 34 corresponds to (3). 

We also treated the cases that p = 2,4 on the same domain 11. Remark that the maximam 
principal guarantees positiveness of all solutions to (7) with even p > 1 except for the trivial 
solution. We derived the estimation results (2) and (4) by computation with bases of N = 140 
and N = 20, respectively. 

In Table 2, one can find a comparison between lower and upper bounds by our method, upper 
bounds by Plum’s formula [18] (Theorem A.2), and the classical bounds by Corollary A.l. 



Figure 1: An approximate solution to (7) in the case that p = 3. 


Table 1: Verification result of the case that p = 3. 


N 

H ( ] -error 

L°°-error 

sup{n (x) 2 \x € fl} 

c 4 (n) 

10 

20 

30 

34 

2.449623e-02 

1.531402e-06 

6.227678e-ll 

2.284208e-12 

3.223795e-02 

1.973787e-06 

8.015445e-ll 

2.939286e-12 

3.458360e-03 

7.071145e-04 

7.070096e-04 

7.070096e-04 

n oo710119440627 
U.Z0 52 444 5 7 6801 o 

0.285244^® 

0.2852444607?]|i 

0.285244460719^ 


Table 2: Upper bounds of the embedding constant by our method, by Corollary A.l, and by 
Theorem A.2. 


c p (n) 

Our method 

Corollary A.l 

Theorem A.2 

c 3 m 

c 4 {n) 

c 5 (n) 

H 0^710766496560 
u.zo/iz 475017617 

0.285244460719!! 

n Q1067136032829 
u ’° iu 58015094169 

0.27991104681667 

0.31830988618379 

0.35780388458050 

0.32964899322075 

0.39894228040144 

0.48909030972535 
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5 Conclusion 


In this paper we proposed a numerical method for deriving a sharp inclusion of the best constant 
in the Sobolev inequality (1). We derived inclusions of the constants by computing the solution 
to problem ( 6 ) with verification. The positiveness of a verified solution to (7) was proved using 
the method proposed in Section 2. The accuracy of all results, e.g., those in Proposition!.. 1, is 
mathematically guaranteed using toolboxes for verified numerical computations [19, 9]. 


A Simple bounds for the embedding constant 


The following theorem gives the best constant in the classical Sobolev inequality. 

Theorem A.l (T. Aubin, 1976 [1] and G. Talenti, 1976 [23]). Letu be a function in H 1 (M n ) (n = 
2,3, • • •). Moreover, let q be a real number such that 1 < q < n, and set p = nq/ (n — q). Then, 


holds for 


— Tp || Vu|| £®(Rn) 



r(i +j)r(n) 

r (f) r ( 1 + ™-?) 


n 


(9) 


with the Gamma function T. 


The following corollary, which comes from Theorem A.l, gives simple bounds for the embed¬ 
ding constant from Hq (O) to L P (Q) for a bounded domain O. 

Corollary A.l. Let O C M n (n = 2,3, - • •) be a bounded domain. Let p be a real number such 
that p € (n/(n — l),2n/(n — 2)) if n > 3 and p € (n/(n — l),oo) if n = 2. Moreover, set 
q = np/(n + p). Then, (1) holds for 

c p (n) = |p| 2 -wt p , 

where T p is the constant in (9). 

Proof. By zero extension outside P, we may regard u € Hq (tl) as a element u € Hq (M n ). 
Therefore, from Theorem A.l, 


IMlLP(fi) < T P ll Vn ll 

Holder’s inequality gives 

ll v «llL ( n) - ( [ \V u ( x )\ q '« dx 


2-q 

2 \ 2 
|1| 2 -« dx 


2-q 

= |H| — 


\Vu (x)|" dx 


that is, 


. 2 -i 


l|V«|| W (n) < l n l 29 II v ' u IIl 2 (Q) • 

where |H| is the measure of H. From (10) and (11), it follows that 

IMIz,p(fi) — 1^1 2q Tp ||Vu|| i 2 (Q) • 


( 10 ) 


( 11 ) 


□ 


6 



Using the following theorem, one can also obtain an upper bound of the embedding constant 
when the minimal point of the spectrum of —A on L/q(D) is concretely computed or estimated. 


Theorem A.2 (M. Plum, 2008 [18]). Let p £ [0, oo) denote the minimal point of the spectrum 
of —A on Hq(LL) for a domain fl C M n . 

a) Let n = 2 and p £ [2,oo). With the largest integer v satisfying u < p/2, (1) holds for 


C'p(fi) 



1 | 2^—3 

2 v r p 

.2 



- — v + 2 
2 


_1 

P 


where i(i- 1 y..(i- I/ + 2 ) = i,f i , = i. 

b) Let n > 3 and p £ [2, 2n/(n — 2)]. With s := n(p _1 — 2 _1 + n -1 ) £ [0,1], (1) holds for 


c p (n) 


( n — 1 \ s _s 

V^(n-2)J P 2 
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